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Abstract 

We show that a method presented in [S.L. Trubatch and A. Franco, Canonical Pro- 
cedures for Population Dynamics, J. Theor. Biol. 48 (1974), 299-324] and later in [G.H. 
Paine, The development of Lagrangians for biological models, Bull. Math. Biol. 44 (1982) 
749-760] for finding Lagrangians of classic models in biology, is actually based on finding 
the Jacobi Last Multiplier of such models. Using known properties of Jacobi Last Multi- 
plier we show how to obtain linear Lagrangians of those first-order systems and nonlinear 
Lagrangian of the corresponding single second-order equations that can be derived from 
them, even in the case where those authors failed such as the host-parasite model. 

PACS: 02.30.Hq, 02.30.Xx, 45.20.Jj, 87.23.Cc 
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1 Introduction 

About 36 years ago, Trubatch and Franco published a paper [20] in which they presented 
an explicit algorithm for constructing Lagrangians of some biological systems, namely the 
classical Volterra-Lotka's model [22], the Gompertz's model [I], the Verhulst's model [21], and 
an host-parasite model [13]. 

Their method for finding a Lagrangian of a second-order equation is, as they state, that 
by Havas [3] who based his method on Helmholtz's work [JJ. Neither Helmholtz nor Havas 
ever acknowledged the use of the Jacobi Last Multiplier in order to find Lagrangians of a 
second-order equation [9^1 [23]. 

Indeed, the method by Trubatch and Franco is based on finding a function / that satisfies 
their equation (6) and is nothing else than the Jacobi Last Multiplier. Because they did not 
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1 An English translation is now available |10] , 
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know the properties of the Jacobi Last Multiplier they were unable to find a Lagrangian for 
the host-parasite model. In fact they found just a linear Lagrangian of this model and stated 
explicitly "In general, there is no relation between the linear Lagrangians of this section and 
the non-linear ones of the previous section for the same model systems." In this paper we prove 
that they were wrong. 

It is interesting to note that the method by Trubatch and Franco for finding linear La- 
grangians is that introduced by Kerner [12]. Again Trubatch and Franco did not realize that 
their key-function W that satisfies their equation (50a) is nothing else than the Jacobi Last 
Multiplied. 

Eight years later, Paine [18] published a paper on the same subject and based his work on 
the method introduced by Kerner [12] . and cited Helmholtz's work [4] as well. Of course, the 
method proposed by Paine is based on a function g that is actually the Jacobi Last Multiplier 
of the two-dimensional systems that he studies. 

Paine posed the following questions: "What are the criteria that a system of ordinary 
differential equations must satisfy to assure the existence of a Lagrangian?" (omissis) "Does 
there exist an algorithm that enables one to construct the Lagrangian from the dynamical 
equations?" 

Strange enough, he did not mention the previous work by Trubatch and Franco [20J. Paine's 
examples are the Volterra-Lotka's model similar to that studied by Trubatch and Franco [20J, 
and two quite trivial linear systems. 

In this paper we show that recognizing that the key-function for finding a Lagrangian is the 
Jacobi Last Multiplier permits to obtain all the results in [20] in a simple and complete way 
and furthermore where Trubatch and Franco fail, namely the model of host-parasite, Jacobi 
Last Multiplier prevails by yielding a suitable Lagrangian. 

The paper is organized in the following way. In the next section we recall the properties 
of the Jacobi Last Multiplier, its well-known connection with Lagrangians of second-order 
equations [23], and the connection with the Lagrangian of systems of two first order equations: 
Noether's theorem [15] is also presented. In section 3 we apply the method of the Jacobi 
Last Multiplier to the same systems as given in [20J, and their equivalent single second-order 
equation. The last section contains some final remarks. 

2 Actually Havas and Kerner never acknowledged each other work, although, at least once, they were pre- 
senting at the same meeting in the same section [2] and [11] , 
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2 The method by Jacobi 



The method of the Jacobi Last Multiplier ( [6] , [7] , [8] , [9] ) provides a means to determine all 
the solutions of the partial differential equation 



Af = ^diixi, 



' Xn) d Xi 







(1) 



(2) 



i=l 

or its equivalent associated Lagrange's system 

¥l _ ?2 __ ¥v 

In fact, if one knows the Jacobi Last Multiplier and all but one of the solutions, namely n — 2 
solutions, then the last solution can be obtained by a quadrature. The Jacobi Last Multiplier 
M is given by 

d(f,ui,u} 2 , ■ ■ .,w„-i) 



d(xi,x 2 , ...,x n ) 



MAf, 



(3) 



where 
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(4) 



and cji, . . . ,u) n —i are n — 1 solutions of (HJ or, equivalently, first integrals of ([2]) independent 
of each other. This means that M is a function of the variables (xi, . . . ,x n ) and depends on 
the chosen n — 1 solutions, in the sense that it varies as they vary. The essential properties of 
the Jacobi Last Multiplier are: 

(a) If one selects a different set of n — 1 independent solutions r/i,... ,r/ n _i of equation ([1]), 

then the corresponding last multiplier N is linked to M by the relationship: 

N = M d( Vl ,...,r )n - 1 ) 
d(ui, . . . ,w n _i) 

(b) Given a non-singular transformation of variables 

T : (xi, X2, • • • , x n ) > (xi, x 2 , . . . , x n ), 
then the last multiplier M' of ^I'i 71 = is given by: 

r d(xi,x 2 , ■■■ ,x n ) 



M = M 



(7 V X 1; x 2' • • • ' x rJ 
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where M obviously comes from the n — 1 solutions of AF = which correspond to those 
chosen for A'F = through the inverse transformation r . 

(c) One can prove that each multiplier M is a solution of the following linear partial differ- 
ential equation: 

i=l 1 

viceversa every solution M of this equation is a Jacobi Last Multiplier. 

(d) If one knows two Jacobi Last Multipliers Mi and M 2 of equation ([1]) , then their ratio is a 
solution u) of (HJ), or, equivalently, a first integral of ([2]). Naturally the ratio may be quite 
trivial, namely a constant. Viceversa the product of a multiplier M\ times any solution 
uj yields another last multiplier M 2 = M\uj. 

There is an obvious corollary to the results of Jacobi mentioned above. In the case that 
there exists a constant multiplier, then any other Jacobi Last Multiplier is a first integral. 

Another property of the Jacobi Last Multiplier is its (almost forgotten) relationship with 
the Lagrangian, L = L(t,x,x), for any second-order equation 

x = <f>(t, x, x) (6) 

i.e. [9| (Lecture 10), [23] 

M = m m 

where M = M(t, x, x) satisfies the following equation 

+ -0. (8) 

Then equation © becomes the Euler-Lagrangian equation: 

d fdL\ dL 

The proof is given by taking the derivative of ([9]) by x and showing that this yields ([8]). If one 
knows a Jacobi Last Multiplier, then L can be obtained by a double integration, i.e.: 



L = J (J MdxJ dx + £i(t,x)x + £ 2 (t,x) 



(10) 



where t\ and £2 are functions of t and x which have to satisfy a single partial differential 
equation related to ([6]) [17J. As it was shown in [17] , ^1,^2 & r e related to the gauge function 
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F = F(t,x). In fact, we may assume 

OX 

h = ^+h(t,x) (11) 

where £3 has to satisfy the mentioned partial differential equation and F is obviously arbitrary. 
In |20| it was shown that a system of two first-order ordinary differential equations 

ui = <f>i(t,ui,u 2 ) 

U2 = <t>2(t,ui,u 2 ) (12) 
always admits a linear Lagrangian of the form 

L = Ui(t,ui,u 2 )ui + U 2 (t,Ui,u 2 )u 2 - V(t,ui,u 2 ). (13) 
The key is a function W such thatH 

^ du2 du\ 

and 

It is obvious that equation (fT5j) is the equation ([5]) of the Jacobi Last Multiplier for system 
(|12|) . Therefore once a Jacobi Last Multiplier M(t, 111,112) has been found, then a Lagrangian 
of system (|12[) can be obtained by two integrations, i.e.: 

L= (J Md«ij u 2 - (J Mdu 2 }iii+g(t,ui,u 2 ) + ^G(t,ui,u 2 ), (16) 

where g(t, ux,u 2 ) satisfies two linear differential equations of first order that can be always inte- 
grated, and G(t, u±,u 2 ) is the arbitrary gauge function^ that should be taken into consideration 
in order to correctly apply Noether's theorem [15] . If a Noether's symmetry 

T = £(i,ui,U2)d t + 77i(i,tii,ti2)d Ul + ri 2 (t,ui,u 2 )d U2 (17) 

exists for the Lagrangian L in (]16p then a first integral of system (I12D is 

dL dL 
-fr- —(771 - fin) - —(7/2 -^2) + G(t,«i,« 2 ). (18) 



We underline that ii\ and «2 always disappear from the expression of the first integral (fT8 
thanks to the linearity of the Lagrangian (I16p and formula (I16p . 



3 In [20] this formula contains an inessential multiplicative constant, namely the integer 2. 
4 The gauge function was not taken into consideration in |20| . 
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3 Some biological examples from [20] 
3.1 Volterra-Lotka's model 

The Volterra-Lotka's model considered in [20] is the following: 

Wl = w\(a + bw 2 ) 
w 2 = w 2 (A + Bw\). 



(19) 



In order to simplify system (|19p we follow [20] and introduce the change of variables 

wi = exp(ri), W2 = exp(r2) (20) 

and then system (fT9j) becomes 



fx = bex.p(r 2 ) + a 
r 2 = B exp(ri) + A. 



(21) 



An obvious Jacobi Last Multiplier of this system is a constant, say 1, and consequently by 
means of (|16|) a linear Lagrangian of system (|21[) is 



t\Ti - r 2 h + 2(— Sexp(ri) + 6exp(r 2 ) - Ar\ + ar 2 ) + — G(t,ri,r 2 ) (22) 



L, 



which (minus the gauge function G) was found in [20J. Moreover we can derive a Jacobi 
Last Multiplier for the Volterra-Lotka system (|19p by using property (b). In fact we have to 
calculate the Jacobian of the transformation (120p between {wi,w 2 ) and (ri,r 2 ) and this yields 
a Jacobi Last Multiplier of system (fT9|) . i.e. 



Mr, 



Mr, 



9(n,r 2 ) 



3(^1,^2) 

Finally, formula f)16|) yields a linear Lagrangian of system ([19 



1 









Wl 






1 









W 2 



W\W 2 



(23) 



Lr^i = log(w 1 )—-log(w 2 ) — + 2(-Alog(w 1 ) + alog(w 2 )-Bw 1 + bw 2 ) + -^-G(t,wi,w 2 ). (24) 
1 J u>2 at 

This Lagrangian was not obtained in |20] , We note that (|19p is autonomous and therefore 
invariant under time translation, namely dt- It is easy to show that the Lagrangian Lr^i in 
(I24p yields a time-invariant first integral through Noether's theorem [15], i.e.: 



9Lr 



+ w 2 - 



9L\u 
dw 2 



A\og(w\) — alog(w 2 ) + Bw\ — bw 2 = const (25) 
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Following [20] we can transform system f)21[) into an equivalent second-order ordinary dif- 
ferential equation by eliminating, say, r\. In fact from the second equation in (|2ip one gets 

n = log {^^) , (26) 
and the equivalent second-order equation in r 2 is the following 

f 2 = -ffeexp(r 2 ) + aj(A-r 2 ). (27) 
A Jacobi Last Multiplier for this equation has to satisfy equation ([5]), i.e.: 

-^(logM) + 6exp(r 2 ) + a = (28) 

namely 

^(logM) + r! = 0, (29) 



by taking into account the first equation in (|2ip . and consequently we get the following Jacobi 
Last Multiplier for equation (|27l) : 

Mi =exp(-n) = — ^— , (30) 
r 2 - A 

the last equality thanks to (|26|) . Then a Lagrangian can be obtained by a double integration 
as in (fT0|) . i.e. 

Li = J B^(r 2 -A)log(^-r2)-r2 + &exp(r2) + ar 2 ) + -^F(t,r 2 ). (31) 

The same Lagrangian (minus the gauge function F) was obtained in [20] . In order to show the 
power of the Jacobi's method we derive at least another Lagrangian for equation (1271) . 

We note that (|27p is autonomous and therefore invariant under time translation. It is easy 
to show that the Lagrangian L\ in (13ip yields a time-invariant first integral, through Noether's 
theorem [15j . i.e.: 

I\ = — ar 2 + r 2 + ^41og(^4 — f 2 ) — 6exp(r 2 ) = const (32) 

As a consequence of the property (d) of the Jacobi last multiplier, the product of a Jacobi last 
multiplier M\ as in (|30p and a first integral I\ as in (|32p of equation (|27p yields another Jacobi 
last multiplier, i.e. 

M 2 = MiJi = -( ar 2 - r 2 - Alog(A - r 2 ) + 6exp(r 2 ) ) (33) 

A — r 2 V / 
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and therefore we can obtain a second Lagrangian of equation (f2T)h i.e. 

L 2 = --(^(Alog(A-r 2 ) -2ar2)(A-r 2 )log(A-r 2 ) 

-{2ar 2 + r 2 )r 2 - 26exp(r 2 ) {[A - r 2 ) \og(A - r 2 ) + r 2 ) 

+b 2 exp(2r 2 ) + 2abr 2 exp(r 2 ) + a 2 r|) + r 2 ) (34) 

This Lagrangian yields another time invariant first integral which is just the square of ii in 
((32]). 

We can keep using property (d) to derive more and more Jacobi last multipliers and therefore 
Lagrangians of equation (j27[) . In fact other Jacobi last multipliers can be obtained by simply 
taking any function of the first integral ii in (|32|) and multiplying it for either M\ in (|3U|) or 
M 2 in ([33j) . and so on ad libitum. 



3.2 Gompertz's model 

The Gompertz's model considered in [20] is the following^: 

wi = wi ( A log ( — - ) + Bw 2 

w 2 = w 2 ^alog (^-^j + bw-^J . (35) 

In order to simplify system (|35p we follow [20] and introduce the change of variables 

tui = mi exp(ri), w 2 = m 2 exp(r 2 ) (36) 
and then system (I35p becomes 

r\ = m 2 B exp(r 2 ) + Ar\ 

r 2 = mi6exp(ri) + ar 2 . (37) 
It is easy to derive a Jacobi Last Multiplier for this system from ([5]), i.e. 

^ log (M [r] ) = -{a + A) => M[ r ] = exp[-(o + A)t] (38) 



and therefore the following Lagrangian 

L[ r ] = exp[— (a + [rir 2 - r 2 fi - 2mib exp(ri) + 2m 2 B exp(r 2 ) + {A- a)rir 2 ] + ^.G(t, n , r 2 ), 

(39) 



d 



°ln [20] some parentheses are missing: an obvious missprint. 
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which (minus the gauge function G) was found in [20j. Then, property (b) yields a Jacobi Last 
Multiplier for the Gompertz's system (|35p . The product of Mr r i in (|38p with the Jacobian of the 
transformation f|36|) between (101,102) and (ri,r2) yields the following Jacobi Last Multiplier 
of system ([35]) . i.e. 



1 J 1 J d(ti;i,t02) 



1 







1 

W 2 



exp[-(a + 



101^2 



(40) 



and consequently a Lagrangian of the original system ([35 
L M = exp[-(a + A)t] 



log(toi)— - log(to 2 )— - 2a log 

102 Wl 



W 2 

m 2 



log(wi) + 2Bw 2 



Wl 



-2bw\ + 2Alog ( — ) log(w 2 ) - (A - a) log(wi) log(w 2 ) 



+ -G(i,«;i,zo2).(41) 



This Lagrangian was not obtained in [20J. 

We can transform system (j37[) into an equivalent second-order ordinary differential equation 
by eliminating, say, r 2 - In fact from the second equation in (j37|) one gets 



r 2 = log 



f\ — Ar\ 
Bm 2 



(42) 



and the equivalent second-order equation in r 2 is the following 



ri 



hn 1 1 exp( ri i + « log ( ) ) ( n - . in ) + .-lr i . ( '13 i 



Using property (b) a Jacobi Last Multiplier for this equation can be obtained. In fact we have 
to calculate the Jacobian of the transformation between (r\,r 2 ) and (ri,ri)> namely (|42p and 
this yields a Jacobi Last Multiplier of equation ([15]) . i.e^ 



M i = M M = exp[-(a + A)t] — ' 



9(ri, ri) rj — Ar\ 

Then a Lagrangian can be obtained by a double integration as in ()10p . i.e. 



(44) 



Li = exp[— (0 + A)i] ^(ri — An) log(ri — An) + mi&exp(ri) — ari log(Bm2) — an J + —F(t, ri ). 

(45) 

The same Lagrangian (minus the gauge function F) was obtained in [20J. 



Of course, we do not consider any multiplicative constants because they are inessential. 
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3.3 Verhulst's model 

The Verhulst's model considered in [20] is the following: 

ib\ = w\(A + Bw\ + f\W2) 

W2 = w 2 (a + bw 2 + h w l)- (46) 
In order to derive a Jacobi Last Multiplier for this system from ([5]), i.e. 

^ log (M H ) + (25 + f 2 )w 1 + (2b + h)w 2 + a + A = (47) 
we assume that Mr w j has the following form: 

M[ w ] = w b ^w b 2 2 exp(6 3 t), (48) 
where 6j, (i = 1,2,3) are constants to be determined. Replacing this My into (|4T|) yields 

, - ^±M 

ABb-Abh + aBb-aBh 
Bb - /i/ 2 

if Bb — /1/2 7^ 0, and therefore if no condition is imposed on the parameters in Verhulst's 
model. Consequently a Lagrangian of system (pIPj) is 



L [w] = exp(M)f^^ +1 -^-^ +1 ^-^ 
\ o\ + 1 o 2 + 



1 &2 + 1 

62+1 h+if^fowi . J^2 , 2a(6 2 + 1) + b 3 \\ d 



that was not obtained in [20] , 

We follow [20] and introduce the change of variables^ 

wi = exp(ri), w 2 = exp(r 2 ) (53) 

and then system (|4"6|) becomes 

fi = A + B exp(ri) + /1 exp(r 2 ) 

r 2 = a + &exp(r 2 ) + / 2 exp(ri). (54) 



7 It is not clear the reason of this change of variables that was performed in [20] 
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We can transform this system into an equivalent second-order ordinary differential equation 
by eliminating, say, r 2 . In fact from the second equation in (|54p one gets 

r 2 = log , (55) 



and the equivalent second-order equation in t\ is the following 

(ah + bri) ri + A 2 b + B exp(2ri)(B6 - /i/ 2 ) - A(a/i + 26ri) 
-exp(n)f/i (aB - f 2 h) + B(2b - h)h - A(2bB - /i/ 2 )U (56) 



1 

h 



Using property (b) a Jacobi Last Multiplier for this equation can be obtained. In fact we 
have to calculate the Jacobian of the transformation between (w\, W2) and (ri,ri), and this 
yields a Jacobi Last Multiplier of equation ([56]) . i.eH 

M 1 = M M = exp(6 1 r 1 + 6 3 t)[ri - A - B exp(n)] fe2 (6 2 + 2){b 2 + 1). (57) 

Then a Lagrangian can be obtained by a double integration as in (jlOp . i.e. 

Li = e X p(6 1 ri)exp(6 3 i)[ri-A-Bexp(ri)] 62+2 + ^ J P(t,ri). (58) 

The same Lagrangian (minus the gauge function F) was obtained in |20j. 
Since a Jacobi Last Multiplier of system (l5"11) is 



M [r] = M [w] 9 ^ llW2) = exp[(6i + l)n + (6a + l)r 2 + M], (59) 
analogously a Lagrangian of system (|54p is 

r 2 ri 



L H = exp((6 1 + l)r 1 + (6 2 + l)r 2 + 6 3 i) 



b x + 1 6 2 + 1 



2 ^ +2 ^ + |^3)) + d G(triir2) (60) 

This Lagrangian (minus the gauge function G) was also obtained in |20j. 



3 Of course, we do not consider any multiplicative constants because they are inessential. 
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3.4 Host-Parasite model 



As stated in [20J, "a simple mathematical model which describes the interaction between a 
host and its parasite and which takes into account the non-linear effects of the host population 
size on the growth rate of the parasite population is given by the equations [13| " 



w\ = (a — bw2)wi 

w 2 = ( A — B— ) W2- 
V w lJ 

As in the previous example it is easy to derive that a Jacobi Last Multiplier is 

exp(Ai) 



Mr, 



2 ' 
W\W^ 



and consequently a Lagrangian of system ([65 
L\ w ] = ex.p[At] 



log(ioi)-^ H 

«j, W1W2 



a B 
2 2 — 



W2 



A 



logOi) 261og(w 2 ) 



w 2 



d 



+ -j^G(t,wi,w 2 ) 



This Lagrangian was obtained in [20J. 

In order to simplify system (|6ip we introduce the change of variable^) 



W2 = T2 exp(At) 



wi = r\ exp(at), 
and then system (I6ip becomes 

h = - 
r 2 = - 



b ex.p(At)rir2 
Bexp(At)r^ 



exp(at)r\ 

Since a Jacobi Last Multiplier of system ([65]) is 

d(wi,w 2 ) 



Mr, 



Mr, 



analogously a Lagrangian of system ([65]) is 



d(ri,r 2 ) rirf 



2 - 



log (n)r 2 h 



nr 2 



2exp(^t) 



br\ log(r 2 ) exp(at) + B d 



(61) 



(62) 



(63) 



(64) 



(65) 



(66) 



+ _G(t,ri,r 2 ). (67) 
r\ exp(at) di 



This change of variables was not performed in |20] 
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This Lagrangian was obviously not obtained in |20| . 

We can transform system (|65p into an equivalent second-order ordinary differential equation 
by eliminating, say, r<i. In fact from the first equation in (|65p one gets 



bexp(At)ri ' 

and the equivalent second-order equation in ri is the following 



n 



bexp(at)ri + B .., 
bexp(at)rf 



r\ + Af\. 



(68) 



(69) 



Using property (b) a Jacobi Last Multiplier for this equation can be obtained. In fact we have 
to calculate the Jacobian of the transformation between (ri,^) and (ri,ri), and this yields a 
Jacobi Last Multiplier of equation ([69jh i.e@ 



Mi = M [? 



<9(ri,r 2 ) _ 6 2 exp(2^)ri 



1 

h 



bexjp(At)rf bexp(At)ri 



bexp(At) 



(70) 



Then a Lagrangian can be obtained by a double integration as in (I10p . i.e. 

L x = bexp(At) log(n) - 6exp(At) log(n) + E^i^l + -^F(t,n). 

exp(atjri ar 



(71) 



This Lagrangian was not obtained in 



4 Final remarks 

As stated by Paine "Among the mathematical results obtained by studying the inverse problem 
of mechanics, it is the explicit algorithms for constructing Lagrangians that offer the model 
builder the most practical benefit." (omissis) "This gives one hope of finding an integral or 
constant of motion for the dynamical system of interest without the hardship of solving the 
system of equations." [18] . 

This paper deals with Jacobi Last Multiplier and its connection to the inverse problem 
of calculus of variation, namely finding one or more Lagrangians for either systems of two 
first-order equations or single second-order equations. 

We could not find better words than those used by Trubatch and Franco at the end of their 
paper "These results clearly show that as the interaction between the populations becomes 
more complex the corresponding Lagrangian becomes more complicated and difficult to find" 
[20] . We hope to address this problem in a future publication. 

10 Of course, we do not consider any multiplicative constants because they are inessential. 
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